AN A POSTERIORI CONDITION ON THE NUMERICAL 
APPROXIMATIONS OF THE NAVIER-STOKES EQUATIONS FOR 
THE EXISTENCE OF A STRONG SOLUTION* 
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Abstract. In their 2006 paper, Chernyshenko et al prove that a sufficiently smooth strong 
solution of the 3d Navier-Stokes equations is robust with respect to small enough changes in initial 
conditions and forcing function. They also show that if a regular enough strong solution exists then 
Galerkin approximations converge to it. They then use these results to conclude that the existence of 
a sufficiently regular strong solution can be verified using sufficiently refined numerical computations. 
In this paper we study the solutions with minimal required regularity to be strong, which are less 
regular than those considered in Chernyshenko et al (2006). We prove a similar robustness result 
and show the validity of the results relating convergent numerical computations and the existence of 
the strong solutions. 
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1. Introduction. For the three-dimensional Navier-Stokes equations the exis- 
tence of strong solutions over an arbitrary time interval is not known. Therefore the 
validity of the results of the numerical solutions of 3d Navier-Stokes equations is not 
obvious. This problem is addressed in Chernyshenko et al (2006) where a rigorous 
relationship between numerical and sufficiently regular exact solutions is given. They 
show for sufficiently smooth initial conditions and forcing functions (data) that al- 
though a priori there is no guarantee of the validity of the numerical solutions, there 
is an a posteriori condition that if satisfied by the numerical results guarantees the 
existence of a strong solution. In this paper we show the validity of the results proved 
in Chernyshenko et al (2006) for the less regular strong solutions which are not cov- 
ered there. 

We will study the Navier-Stokes equations in their functional form. For a bounded 
domain Q we let TL be the space of divergcncc-frcc smooth vector- valued functions on 
Q with compact support and zero average and define 

H = closure of H in [L 2 (n)} 3 , 
V = closure of H in [H 1 ^)] 3 . 

We use the same notation H and V for the similar spaces of periodic functions over 
the periodic domain Q. Then the Navier-Stokes equations in their functional form 
are written as (Constantin and Foias 1988, Robinson 2001) 

dz£ 

— + vAu + B(u, u) = f, with it(0) = u (1.1) 

where Au = — IIA, B(u, u) = H(u.V)u with II the orthogonal projection from L 2 into 
H. We will consider this equation with the following cases for the data 

(a) 'minimal regularity' when u € V and / e L 2 (0, T; H) n L^O, T; V) 

(b) 'second order regularity' when u € V 2 and / e L 2 (0, T; V) n L x (0, T; V 2 ) 7 
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where V m = H m n V. For the periodic case we know that D(A m / 2 ) = H m n V 
for all m and therefore we define the norm on V m as \\u \\ m = \A m / 2 u\. We denote by 
(•, •) and | • | the inner product and norm on H. 

We will show that the strong solution of with the data introduced in (a) or 
(&), if it exists, remains strong if the changes in the initial condition and forcing func- 
tion are small enough. The exact conditions required for these changes to be 'small' 
are given in theorems l4.ll and l4.2l For example in the case of minimally regular strong 
solutions we require 



\D(u -v )\ + / \Df(s)-Dg(s)\ds 
Jo 

, i/4 / k 2 r 

k\2 7T ) exp HrX 



<lb? expU- \Du(s)\ 4 + -\Du(s)\\Au(s)\ds 



v 



We then, in theorems 16.11 and 16.21 use these robustness results to find an a posteriori 
condition that if satisfied by sufficiently refined numerical approximations, implies 
the existence of a strong solution. We also show that if a strong solution exists the 
Galerkin approximations converge to it and then use this to prove that the existence 
of a strong solution can be verified by the Galerkin approximations. In the last section 
we consider a channel flow as a physical example that can be described by the Navier- 
Stokes equations with the conditions introduced above. For this example we will show 
how the results of this paper can be applied to the Galerkin approximations to verify 
the existence of a strong solution. 

The results we prove here for a strong solution with lowest regularity hold in a 
general bounded domain as well as in the absence of boundaries unlike the results for 
more regular solutions which are proved only for the equations in a periodic domain 
or the whole space. 

2. General ODE lemma. We first prove an ODE lemma which will be used in 
dealing with the differential inequalities that appear in the proofs. We consider the 
differential inequality 

^<8{t)+ay n with y(0) = y >0 
at 

and find the conditions on yo, S(t) and a that ensure that y(t) exists on a finite time 
interval [0, T]. This lemma is a generalization of the result obtained for n = 2 in 
Chernyshenko et al (2006). 

Lemma 2.1. Let T > 0, a > and n > 1 be constants and let S(t) be a non- 
negative continuous function on [0,T]. Let y satisfy the differential inequality 

dv 

tf < S(t) + ay n with y(Q) = y > (2.1) 



dt 



and define 



(*) If 



V = Va+ I S(s)ds. 



V< [(n - l)oT]V(«-i) (Z2) 



then y(t) remains bounded on [0, T] 



(ii) y{t) — > uniformly on [0, T] as 77 — > 0. 
Proo/. 

We first consider the following differential inequality 

z < az" with z(0) = rj (2.3) 

and show that 

sup y{T) < sup z(T). (2.4) 

Si s 2 

where Si and S2 are the sets of all possible solutions of inequalities (|2.1[) and (|2.3[) 
respectively. 

Since y, i, y(0) and z(0) are non-negative, the suprema of y(T) and z(T) happen 
when 

y=5(t)+ay n , 
z= az n 

for all t € [0, T]. In this case, for the difference w = y — z we have 

w = 6(t) + awy2( n ~ )y n - 1 - k z k with w(Q) = - S(s) ds. 

Since y(t) and z(t) are greater than zero and assuming they remain finite over [0, T], 
there exists some M > such that 

w < aMw + 5(t) 

and therefore by Gronwall's inequality 

w{t)<w(0)e aMt + f 5(s)ds 
Jo 

= {l-e aAIt ) f S(s) ds. 
Jo 

This implies w(t) < for any t G [0,T] and (|2~4")l follows. 

The inequality (12. 3p can be written as 

z 

— < a 

z n - 

Integrating both sides from to t yields 

Tj 

Z{T) ~ (-aT(rz - 1)(7 7 )— 1 + l)V«-i " (2 ' 5) 
Since y(t) < y(T) < z(T), y(t) remains finite on [0, T] provided that 

aTtn-l)^)"- 1 < 1 



which yields (|2?2j) . 

Furthermore, it is clear from (|2.5|) that z(T) — > as 77 — > 0, from which it follows 
that y(r;) — > uniformly on [0, T}. □ 
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3. The inequalities for the nonlinear term. For the nonlinear operator, we 
will use the following inequalities in this paper 

\(B(u,v),Aw)\< k\Du\ \Dv\ 1 / 2 \Av\ 1 ' 2 \Aw\ (3.1) 
\(B(w,u),A 2 w)\<c\\u\\ 3 \\w\\ 2 (3.2) 
\(B(u,w),A 2 w)\<c'\\u\\ 3 \\w\\l (3.3) 

The first inequality holds for both Dirichlet and periodic boundary conditions (Con- 
stantin and Foias 1988) while the other two are valid only in the absence of boundaries 
with c and d independent of the size of the domain (Kato 1972). The constant k for 
the periodic and no-slip domain and also the whole of M 3 , is independent of the do- 
main f2. In a general domain with nonzero boundary conditions however, it depends 
on the regularity properties of the boundary of Q. To see this we use the Holder 
inequality to write 

3 

\B(u,V,Aw)\< ^ hi\\L4 D i V j\\L^\AWj\. 

By the Sobolev inequality 

HUe<c|D«| (3.4) 

we have 

II^IUs < WDvW^lDv] 1 / 2 < c^Av^Dv^ 2 . 
Therefore we can write 

3 

\B(u,v,Aw)\ < Y W^\\l4 d i v j\\lM w j\ 

<9\\u\\ L e\\Dv\\ L 3\Aw\ 

< 9c 3 J 2 \Du\ \Dv\ 1 ' 2 \Av\ 1 / 2 \Aw\ 

and so k = 9c^ 2 . For no-slip boundary conditions and the whole of R 3 , c s , the 
constant of the Sobolev inequality (|3.4p . does not depend on f2 (Ziemer 1989). For 
a bounded domain with non-zero boundary conditions, c s depends on the regularity 
properties of d£l but is independent of the size of (Adams 1975). Therefore for the 
cubic domain of the periodic boundary conditions also, it does not depend on the size 
of the domain. In fact, the cubic domain has the strong Lipschitz property and for 
such a domain Adams (1975, Lemma 5.10) has shown that c s — Ay/2. Similarly, c and 
c' depend on the constant of the Sobolev inequality ||?/|| i 6/(3-2(s) < c Sj fc||it||#jc which 
is again independent of the size of the domain (Adams 1975 and Ziemer 1989). 

In obtaining the robustness conditions in the next section, we keep track of the 
constants fc, c and c' when they appear, so that the value of the constant coefficients 
in the robustness conditions could be computed if desired. 

The inequalities (|3.1|) - (|3.3|) are not as elegant as the inequalities available for 
the more regular solutions considered in Chernyshenko et al (2006) and this is why 
different proofs are needed in the less regular cases we are studying here. 

We note here that assuming a strong solution u e L°°(0,T;V) n L 2 (0,T;V 2 ) 
exists, from (|3. 1[) we have B(u,u) S L 2 (0,T; H) and therefore (jl.ip implies that 
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du/dt G L 2 (0,T;H). Having u G T^ 2 and / G I?(0,T;V) in case (b), from the 
regularity result for periodic domains proved in Constantin and Foias (1988) we know 
that the strong solution u is in fact more regular and an element of L°°(0,T; V 2 ) n 
L 2 (0,T;F 3 ). 

4. Robustness of strong solutions. Here we show that as long as a strong 
solution exists for some specific initial data and forcing function, the equations with 
sufficiently close data also have a strong solution. A similar result about the robust- 
ness of strong solutions with respect to changes in the forcing function is proved by 
Fursikov in his 1980 paper. For the three-dimensional Navier-Stokes equations with 
initial condition u G V m where m > 1/2, he shows that the set of forcing functions 
for which a strong solution exists is open in L q (0,T; V" 1 ^ 1 ) with q>2. 

For both minimally regular solutions (a) and the more regular case (b) the condi- 
tion we obtain here for the data depends explicitly on the viscosity coefficient unlike 
the more regular cases considered in Chernyshenko et al (2006). The robustness result 
we prove for the minimally regular strong solutions holds in a general bounded do- 
main. For the second order regular solutions however, we need to restrict the domain 
to be periodic or the whole of R 3 since the inequalities (|3 . 2|) and (|3.3[) only hold in a 
periodic domain (Constantin and Foias 1988) or the whole of R 3 (Kato 1972). 

(a) Strong solutions with minimal regularity. To prove the robustness with 
respect to the data we write the governing equation for the difference between the 
solution of the equations with nearby data and the strong solution and find a bound 
on the norm of this difference. 

Theorem 4.1. Let f G L 1 (0, T, V) fl £ 2 (0, T; H), u eV and u e L°°(0,T;V)n 
L 2 (0,T;V 2 ) be a strong solution of 

du 

h Au + B(u, u) — f with u(0) = uq. 

dt 

If 

|£>(«o-«d)|+ / \Df(s)-Dg(s)\ds 
Jo 

1 / ;/ 3 \ 1/4 / k 2 r T 27k 2 1 1 \ 

<- k {^) -p(- Y / o ~\Du^ + -\I^(s)\\Au(s)\^ (4.1) 

then the solution of 

dv 

h vAv + B(v, v) — g with v(0) = vq 

dt 

is also a strong solution on [0, T] with the same regularity as u. 

Proof. By the local existence of strong solutions (Constantin and Foias 1988, 
Temam 1995) we know that there exists T* > such that v G L°°(0,T';V) H 
L 2 (0,T';V 2 ) for any T 1 < T* . We consider T* the maximal time of existence of 
the strong solution v meaning that limsup t ^ T » \Dv\ = oo. In the following argument 
we assume T* <T and deduce a contradiction. 

We consider w = u — v which satisfies 

h vAw + B(u, w) + B(w, u) — B(w. w) = f — g. 

dt 
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Over the time interval t € [0,T') for any T < T* we have dv/dt e L 2 (0,T';H) and 
since T* < T, du/dt € L 2 (0,T'; H). Therefore taking the inner product of the above 
equation with Aw and using (|3.ip we can write 

~\Dw\ 2 + v\Aw\ 2 < k\Du\\Dw\ 1/2 \Aw\ 3 ^ 2 + k\Dw\\Du\^ 2 \Au\ 1/2 \Aw\ 

+k\Dw\ 3 / 2 \Aw\ 3/2 + \D(f - g)\ \Dw\. 

We then use Young's inequality to remove \Aw\ (which causes the appearence of v in 
the coefficients) and obtain 

i a /97k 4 h 2 \ 97k 4 

n^\DM 2 < -TTlDuf + -\Du\\Au\ \Dw\ 2 + —^iDw] 6 + \D(f - g)\ \Dw\. 



2dt' 1 - V 4^ 3 1 1 2v l 11 7 1 1 4z/ 3 

Dividing by \Dw\ (assuming \Dw\ is nonzero. If not, replacing \Dw\ by \Dw\ + eo 
with eo > will lead us to the same result after tending eo to zero at the end) we get 

H /97k 4 k 2 \ 97k 4 

±\Dw\ < {-£r\ Du \ + ^\Du\\Au\) \Dw\ + -^-flDwf + \D(f g)\. 

Now letting 

and setting y(t) = \Dw{t)\e~^^ the above inequality can be written as 

^<ay 5 + S(t), 

where a = ^e 4f3{T) and 8{t) = \Df{t) - Dg{t)\. 

By lemma [27X1 if the condition (|4.1[) is satisfied, |Dio(f)| is uniformly bounded on 
[0, T*). This implies that \Dv(T*)\ < oo which contradicts the maximality of T*. It 
follows that v(t) is a strong solution on [0, T]. □ 

(b) Strong solutions with second order regularity. We follow the same 
approach as in the minimally regular case to show the robustness for a strong solution 
with second order regularity. Here also the inequalities we need to use for the nonlinear 
operator to find a bound on the V 2 -norm of the difference between nearby solutions 
depend on the l^ 3 -norm of this difference and therefore the robustness condition on 
the data depends explicitly on the viscosity coefficient. However, its dependence is 
weaker in this case. 

Theorem 4.2. Let f e L x (0, T, V 2 )C)L 2 (0, T; V), u Q e V 2 and u e L°°(0, T; V 2 )D 
L 2 (0,T;V 3 ) be a strong solution of 

du 

— + Au + B(u, u) = f with u(0) = u . 



If 



£\A(f-g)\dt < I^exp^ 



\A(u -v a )\+ I \A(f-g)\dt < -W^expl- Ac + c')\\u\\ 3 dt ) , (4.2) 



then the solution of 



— — h Av + B(v, v) = g with v(0) = vq 

is also a strong solution on [0, T] with the same regularity as u. 
Proof. As before, the difference w — u — v satisfies 

— — + vAw + B(u, w) + B(w, u) + B(w, w) = f — g. 

Taking the inner product of the above equation with A 2 w and using (|3.2p and (|3.3p 
we obtain 

~\\wh + u\\w\\i<(c + c')\\uf^ 

We apply Young's inequality to the second term on the right hand side and get 

^IHIl<(c + c')lH|3|Hll + ^||Hlf + ll/-flll2|IHl2. 

Dividing the above inequality by \\w\\2 yields 

^\\wh < (c+0||«|| 3 Hh + ^IHII + II/- -9 1 1 2 - 
Now we consider 

y(t) = \\whe-M 

with P(t) = J (c + c')||w||3 dt. The equation with this new variable becomes 

where a = j^e 213 ^ and 5 = \\f — <7 1 1 2 . Using lemma EOl ||io||2 is bounded if the 
condition (|4.2[) is satisfied. □ 

5. Convergence of Galerkin approximations. Here we show that if a strong 
solution with minimal or second order regularity exists, Galerkin approximations con- 
verge to it. Similar results about convergence of various numerical method assuming 
the existence of a strong solution, are given for finite element methods by Heywood 
and Rannacher (1982), for a Fourier collocation method by E (1993) and for a nonlin- 
ear Galerkin method by Devulder, Marion and Titi (1993). Here as in Chernyshenko 
et al (2006), we make no assumption on the regularity of the Galerkin approximations. 

For the minimally regular case the result we prove here holds for the solution of 
the equations over a general bounded domain. For the second order regular solutions 
we need to use inequality (|3.3[) which is valid only in a periodic domain or the whole 
of R 3 . 

In the following theorems we let P n be the orthogonal projection in H onto the 
space spanned by the first n eigenf unctions of the Stokes operator A, {u>j}™ =1 , or- 
dered so that their corresponding eigenvalues satisfy < Ai < A2 < . . . . Since the 
eigenfunctions of A are smooth (Constantin and Foias 1988), for any u € V m , m > 
(with V° = H) we have u n = P n u = Y^j=i{ u i w i) w 3 e V m ■ 

We note that what we obtain here about the convergence of Galerkin approxima- 
tions would be useful even if the existence of regular solutions was known. 
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(a) Minimal regularity. The following theorem, like the robustness theorem in 
the minimal regularity case, holds in sufficiently smooth bounded domains as well as 
in the absence of boundaries. 

Theorem 5.1. Let u a e V, f e L 2 {0,T;H) and u e L°°(0, T; V) n L 2 (0, T; V 2 ) 
be a strong solution of the Navier-Stokes equations 



du 
dt 



+ vAu + B(u, u) = f{t) with tt(0) = uq. 



Then u n , the solution of Galerkin system 



du 



+ vAu n + P n B{u n ,u n ) = P n f{t) with u n (0) = P n u , 

converges strongly to u in both L°°(0, T; V) and L 2 (0, T; V 2 ) as n — > oo. 
Proof. We consider w n = u — u n which satisfies 



(5.1) 



dw n 
dt 



vAw n + P n Biu,w n ) + P n B(w ni u) - P n B(w n ,w n ) = Q n f{t) - Q n B(u,u). 



Letting q n = u — P n u and taking the inner product of the above equation with Aw n 
while noting that {P n B(u, v), Aw n ) = b(u,v, P n Aw n ) we obtain 

~\Dw n \ 2 + »\Aw n \ 2 < k (\Du\ IDw^Aw^ 3 ' 2 

+ \Dw n \ \Du\^ 2 \Au\^ 2 \Aw n \ + \Dw n \ 3 l 2 \Aw n \ 3 ' 2 
+\Du\ 3 / 2 \Au\ 1 / 2 \Aq n \) +\f\\Aq n \. 

Applying Young's inequality and redefining the constant k gives 



dt 



\Dw n \ 2 <k 



1 



1 



1 



Duf + -\Du\ \Au\ \Dw n \ 2 + -^\Dw n 



Now letting 



and 



+\Du\ 3 ' 2 \Au\ 1 / 2 \Aq n \} + \f\\Aq n \. 
0(t) = kj ^|£>u| 4 + \\Du\ \Au^j ds 



y n (t)=e-^\Dw n \ 2 
the above inequality can be written as 

^<ayl + S n (t), 

with a = £e 2 « T ) and 6 n (t) = \f\ \Aq n \ + klDu^Auf^Aq* 
Using the Holder inequality we have 



Jo 



S n (s) ds < 



t \ 1/2 

2aA , f. 



\f\ 2 ds 



T \ V2- 

\Du\ 3 \Au\ds 



1/2 



\Aq n \ 2 ds 



Since u e L 2 (0 : T;V 2 ), u(s) e V 2 for almost every s G [0,T] and therefore <z„(s) = 
Q„m(s) — > in V 2 as n — > oo and |>l<7„(s)| 2 < |Au(s)| 2 for a.e. s e [Q,T]. Hence the 
Lebesgue dominated convergence theorem for each s implies that J \Aq n \ 2 ds tends 
to zero as n — > oo and therefore 

T 

5 n (s) ds — > 

as n — > oo. Since y„(0) — ► as n — * oo, lemma [2TT1 shows that y„, and hence 1-DuVil 2 , 

converges to zero uniformly on [0, T] as n — ► oo. 

□ 

(b) Second order regularity. Here also, in a similar way to the previous case, 
we try to find a bound on the l/ 2 -norm of the difference between a second order 
regular strong solution and Galerkin approximations. However here the more regular 
solution space makes the proof easier. 

Theorem 5.2. Letu G V 2 , f G £ 2 (0, T; F)nL 1 (0, T; F 2 ) ondue L°°(0, T; F 2 )n 
L 2 (0,T; V 3 ) 6e a strong solution of the Navier-Stokes equations 

du 

— + vAu + B{u, u) = f(t) with u(0) = Uq. 

Then u n , the solution of Galerkin system 
du 

—— + vAu n + P n B(u n ,u n ) = P n f{t) with u„(0) = P n u , (5.2) 
dt 

converges strongly to u in both L°°(0, T; V 2 ) and £ 2 (0, T; V 3 ) as n — ► oo. 
Proof. Let w n = u — u n . Then w n satisfies 

— — + vAw n + P„B(u, w n ) + P n B(w n ,u) + P n B( 

w m W n ) — Qnf QnB(u, It) 

dt 

Taking the inner product of the above equation with A 2 w n and using (|3.2p and 
we obtain 



Id 

2dt 

+ \\Qnf ~ Q n B(u,u)\\ 2 \\w n \\2 



77 — \\ w n\\l + HWIi < (c + c / )||tt|| 3 ||ix; n ||2 + c||ty n || 3 ||«; n ||l 



Now we use Young's inequality to remove the dependence on 1 1 it?n 1 1 3 and then divide 
by || w || 2 to obtain 

d c 2 

TrlWh < (c + c' ) 1 1 -ix 1 1 3 1 1 zi7„ 1 1 2 + -HI^Hl + WQnf - Q n B(u,u)\\ 2 , 
dt \v 

in which the coefficient of ||w n ||2 does not depend on n. 
Letting 

y n = \\w n he-M, 
with f3(t) — exp (^f*(c + c')\\u(s)\\ 3 ds^j yields 

Vn < 5(t) + ayl, 



where a = f I7 e 2 ' 3(T) and 8{t) = \\Q n f - Q n B(u, u)\\ 2 . So by lemma[2j]if 

V = \\Q n uo\\m + / \\Qnf - QnB(u,u)\\ 2 ds -> as n -> oo, 
Jo 

then y n (t) — > uniformly on [0, T] as n — > oo. 
By (I5.2[) ||Q n Uo||m — > as n — > oo. We know that 

||B(u,«)||2<c||u|| 2 ||u||3 (5.3) 

(Kato 1972, Constantin and Foias 1988) and therefore f(s) - B(u(s),u(s)) G V 2 
for almost every s G [0, T]. So since {u>j}™ =1 form a basis for y 2 as well as H, 
\\Q n (f(s) - B(u(s), u(s))) ||a -> and 

||Q„ (/(a) - «(«))) || a < ||/(s) - B{u{s), u(s))\\ 2 

for almost every s € [0, T] . Therefore by the Lebesgue dominated convergence theorem 
it follows that 

f ||Q„(/(«)-S(u(s),u(«)))|| 2 -0 
Jo 

and the result follows. 
□ 

6. Numerical verification of the existence of a strong solution. Here we 
show that the existence of minimal and second order regular strong solutions can be 
verified via computations using sufficiently refined Galerkin approximations. 

Theorem 6.1. i) Consider the Navier- Stokes equations 

dzt 

\- vAu + B(u, u) = f with u(0) — uq (6-1) 

dt 

with Mo G V, f G L 2 (0,T; H) n L 1 (0,T; V), that hold in a bounded domain ft with 
sufficiently smooth boundary or in a periodic domain. 

Let v G L°°(0, T; V) n L 2 {Q, T\ V 2 ) be a numerical approximation of u satisfying 

^ + vAv + B(v, v) G L x (0, T; V) n L 2 (0, T; fl") 
dt 

and 

\Dv{Q) - Du \ + [ T + vAv{s) + B(v(s),v(s)) - ds 

Jo ds 

1 / u 3 \ 1/4 / k 2 r T 27k 2 1 1 \ 

<- k {^) ^[-yI — ^W-)| 4 + -W-)ll^(.)|d S J.(6^) 

T/ien i/ie Navier-Stokes equations \6.1\) have a strong solution u G L°°(0,T;V^) PI 
L 2 (0,T;V 2 ). 

ii) Let u be a strong solution of i6.1\) . Then there exists an N such that the 
Galerkin approximations u n satisfy the inequality I6'.£l) for all n > N . Therefore in 
view of part (i), u passes the a posteriori test as a solution approximated by u n i.e. the 
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existence of a strong solution with minimal regularity can be verified by the Galerkin 
approximations. 



Proof. 

i) Considering 



v is a strong solution of 



dv 

g = — + vAv + B{v,v), 



dv _ 

— — h vAv + B(v, v) = q. 

dt v ' 

So by theorem 14.11 if the inequality (|6.2[) holds, the solution of the equations with 
nearby data (u Q , f) is a strong solution. 

ii) The strong convergence of u n to u in L°°(0,T;V) and L 2 (Q,T;V 2 ) is guar- 
anteed by theorem 15.11 Therefore the right hand side of (|6.2[) is bounded below for 
every n. It follows from (15. 1|) that |Du„(0) — Duo\ — > as n — > and 

+ i/Ait„ + B(u n ,u n ) - f = Q n [B(u n ,u n ) - /]. 

Therefore it only remains to show \\Q n [B(u n (s),u n (s)) — f]\\\ converges to zero as 
n — > oo for a.e. s e [0, T], it is then clear that (|6.2| will be satisfied for all n sufficiently 
large. 

For the nonlinear operator we have 

||B(u„,u n )||i = \J {D{u n ■ V)u n , D{u n ■ Vjun) < c (\\Du n \\ 4 L4{n) + \\u n \\ 2 L ^ {n) \Au n \ 2 ^j 
Using the Sobolev inequaities we can write 

\\Du n \\ L i {n) < c\\Du n \\ L 6^ < c\Au n \ 

and 

||Mn||i~ < c|Am„|. 
Therefore ||B(u n ,zt n )||i < c|Au„| 2 which implies 

f(s) - %(«), u n (s)) G V for a.e. s £ [0, T] 

and we will have 

\\Q n {B(u n {s),u n (s)) — /(s))||i — > as n — > oo and for a.e. s G [0, T] 

and 

||Q n (BK(s),n„(s))-/(s))|| 1 < \\B(u n (s),u n (s))-f(s)\\ for a.e. sG [0,T]. 
By the Lebesgue dominated convergence theorem we conclude 

T 

\\Q n (B(un(s), u n (s)) - /(s))||ids — > as n > oo. 
11 
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□ 

A similar result holds for the strong solutions with second order regularity in a 
periodic domain. 

Theorem 6.2. i) Consider the Navier- Stokes equations 

du 

— + vAu + B(u, u) — f with u(Q) — uq (6-3) 

with uq € V 2 , f G L 2 (0,T; V 1 ) fl L 1 (0,T; V 2 ), that hold in a periodic domain. 

Let v € L°°(0, T; V 2 ) n L 2 (0, T; V 3 ) be a numerical approximation of u satisfying 

-^- + vAv + B(v, v) S L x (0, T; V 2 ) n L 2 (0, T; 7) 
di 

and 

ll«(0) - u4l + £ II ^(s) + vAv(s) + B(v(s),v(s)) - f(s)\\ 2 d s 

< i^cxp (-J\ c + c')\\u n \\ 3 dtj . (6.4) 

Then the Navier-Stokes equations H6.3\) have a strong solution u G L°°(0, T;V 2 ) fl 
i 2 (0,T;T/ 3 ). 

ii) Lei u be a strong solution of i6.S\) . Then there exists an N such that the 
Galerkin approximations u n satisfy the inequality |6'.4| i /or all n > N . Therefore in 
view of part (i), u passes a posteriori test as a solution approximated by u n which 
means that the existence of the strong solution u can be verified using the Galerkin 
approximations. 

Using the robustness result for the strong solutions with second order regularity, 
theorem 14.21 the proof of the above theorem closely parallels that of theorem 16.11 

7. A physical application: channel flow. We considered the Navier-Stokes 
equations for a flow in a domain with non-moving boundary and with the forcing 
function in the form of a body force. A physical situation for which the results of 
this paper and in particular the a posteriori test of theorem 16.11 can be applied is 
a channel flow in the domain < x < L Xl < y < 1, < z < L z with the 
velocity field u(x,y, z,t) considered to be periodic in x and z with periods L x and 
L z respectively and zero at y — and y = 1. The body force g is assumed constant 
and equal to (1,1,1). Here we show how to write the inequality (16. 2\ for the Galerkin 
approximations, u n , of the solution u of this problem. 

By theorem 16.11 we know that if for some n the Galerkin approximation, u n , 
satisfies 

f T du (s) 

\Du n (0) - Du Q \ + I \\—^ J -+uAu n {8) + B{u n {8),u n {a))-f(a)\\ 1 d8 

Du n (s)\ 4 + -| Du„ (s)| \Au„ (s)\ds) 




(—) 

k \27TJ 



then a strong solution u exists and the Galerkin approximations converge to it. To 
check this inequality for a Galerkin approximation of the above channel flow example 
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we need to compute Au ni B(u n , u n ) and / = Tig. 
The functions 



2-rri(-p^x + !p-z) ■ i r \ 



with k = (fci, fc2 ; ^3)1 form an orthogonal basis for the space of L 2 -functions on the 
cubic domain introduced above, which have periodic bounday values in x and z di- 
rections and are zero on y = and y = 1. We can therefore define the Galerkin 
approximation u n in terms of u>k as 

n 

U n = ak^k, with V • U n = (7-1) 

k=n 

where a k = (a lk (t),a 2k (t),a 3 ^(t)), a-(_ felife2) _ fc3 ) = o^ kl>k2tk ^, n = (-rc,0,-ra) and 
n = (n, n, n). Since u n is divergence- free we have 

n k k 

^ 27ri(— aik + — «3k) sm{irk 2 y) + nk 2 a 2k cos(irk 2 y) = 0. 

k=n 

After expanding cos(jrk 2 y) in terms of sine functions we obtain 

[ fc 2+" j 

^ ( | aik + | a3k)+ ^ ^ + 1 - fc2)( 2fc2 - 1 2^i + 27TT )Q2k = - 

Defining 

[ fc 2+" j 

and having (|7.1[) we obtain 

1.2 1.2 1,2 



and 



k=n x z 



2„ k+ L^^J 



B(u„,u„)=2^ ^ ^ <j (a km -j • j) ( Oj - 7^(«j ' k) 

k=2n j=k_ m=L ^2_tij I V l k l / 

1/1 1 1 1 



T:\2m~ 2j 2 + l 2m~2j 2 -2k 2 + l 2m + 1 2m - 2k 2 + 1 



where k m = [k\,2m + 1 — k 2l k 3 ), k = (min{fci, 0}, min{fc2, 0}, min{fc3, 0}) and 
k + = (max{fci, 0}, max{/c2, 0}, max{fc3, 0}). 

Now we note that / = (1,0,1), since / — g = (0,1,0) is perpendicular (with 
respect to L 2 -norm) to all functions in H. Therefore using the expansion of u n and 
Au n with respect to {wk}jj =n and B(u ni u n ) with respect to {wk} k ™ 2no we can wr ite 
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(noting that the coefficients of % when ki, k 2 or /c 3 are less than — n or bigger than 
n are zero for u n and Au n ) 

d 2n 
E n = —r^ + vAu n + B{u n ,u n ) - f= ^2 /?kw k - (1,0,1) 

k=2n n 

and therefore 

' '' " ' k=2n 



and 

2n 

1 4- uAii... 4- /■': ,/ /, : - fll? = 2r,_7,. 

di 



^ + j/Au,, + %, u n ) - /||? = 2L X L, £ (/3 2 k + /3| k + /3 3 2 k )(fc! + k 2 + fc 3 ) 2 . 



k=2n 

The norms of Du n and Aw n are computed as 



|Du„| 2 = 2L X L Z («ik + a 2k + «!k)(ki + ^2 + fc 3 ) 2 , 

k=n 

|Au„| 2 = 2n"L x L z J2 « + < + «L)(|| + f + %?■ 



k=n 



The remaining term on the left-hand side of (|6.2p , |Du n (0) — DuqI, is obtained in a 
similar way to |Z?u„|. 

Finally we need the value of the constant k in the inequality (|6.2[) . The cubic do- 
main here also has the strong local Lipschitz property and as it was shown in section 
3, for such a domain k = 72(2 3 ' 4 ). 



8. Conclusion. We extended the results of Chernyshenko et al (2006) to the less 
regular strong solutions not considered in their paper. Unlike their results, the ro- 
bustness conditions on the data that we obtained for the minimally and second order 
regular strong solutions in theorems 14.11 and 14.21 depended explicitly on the viscosity 
coefficient in a way that as the viscosity coefficient became smaller, the conditions 
(|4.ip and (|4.2p became more restrictive which obviously is not desirable. Moreover, 
since this is due to the less elegant inequalities available for the trilinear form in the 
less regular spaces considered here, it does not seem that we can remove the viscosity 
coefficient from the robustness conditions. 

On the other hand, the results proved here for the minimally regular strong so- 
lutions are valid in a sufficiently smooth bounded domain as well as in the absence 
of boundaries while in the case of the second order regular solutions considered also 
here and the more regular ones studied in the above mentioned paper the results are 
restricted to a periodic domain or the whole of R 3 . The reason is that the inequalities 
we need to use in these cases for \(B(u,v),v) m \ when m > 2 are valid in the absence 
of boundaries where we know that swapping the action of the Leray projector and the 
Laplacian operator is possible. 

It would be interesting to see if there exists a big enough m for which it is possi- 
ble to swap the action of the Laplacian operator and the Leray projector in D(A m ^ 2 ) 
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over a bounded domain. If such a finite m exists we can extend the robustness results 
for regular enough solutions to the more physically relevant bounded domains. Fur- 
thermore, for two dimensional bounded domains equality of the Stokes and Laplacian 
operators in D{A m / 2 ) can be used in obtaining a better bound on the attractor di- 
mension in L 2 . 
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